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Abstract 



Since the electromagnetic current for a pion system coincides with the third 



'"qJ component of the isovector current, the isospin angular oscillation of a co- 

Qh! herent field can be a significant source for the electromagnetic emissions. We 



study the characteristic dilepton and photon emissions from the classical pion 



field oscillation in the QCD vacuum. The general analytical solution obtained 
in the nonlinear a model is used to calculate the electromagnetic current den- 
sity, which exhibits a light-front singularity and decreases rapidly as inverse 
square of the proper time due to a longitudinal expansion. The momentum 
and invariant mass spectra of the direct photon and dilepton are found to be 
a sensitive probe of the space-time evolution of the chiral condensate field. 
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Recently, there has been some interest in the formation of a disoriented chiral conden- 
sate (DCC) in high energy hadronic and heavy-ion collisions (]]]. Preliminary studies have 
suggested that there may be a coherent pion field oscillation following a nonequilibrium 
second order phase transition which may lead to characteristic soft pion production 0|§. 
According to a "Baked Alaska" scenario suggested by Bjorken, Kowalski and Taylor |TJ], the 
collision debris in a high energy central event form a "hot" source on the light front and 
move outward while emitting quasi-goldstone bosons in the cool interior of light cone. These 
pions are not assumed to reach a thermal equilibrium, rather, they behave collectively and a 
coherent field description is more relevant. The space-time evolution of the condensate field 
is mainly determined by the classical equation of motion for a low energy effective theory 
such as a linear or nonlinear a model. Eventually, the "disoriented vacuum" relaxes back to 
normal vacuum by emitting physical pions. It is, however, difficult to obtain the information 
on its space-time evolution from the final hadron spectrum. 

In contrast, the electromagnetic signature is known to be an ideal probe of a dense 
hadronic matter jJJ owing to the fact that it escapes the strong interaction region once 
produced without further final state interactions and thus carries the information on the 
early dynamical evolution. In this Letter, we shall study the dilepton and direct photon 
production from the classical pion field in the context of a disoriented chiral condensate 
(DCC) or most generally a nonequilibrium pion cloud. We shall develop a general formalism 
for the dilepton and photon emissions in the presence of a classical electromagnetic current. 
As an example, we calculate the dilepton and photon spectra in an analytical solvable model 
where a general class of solutions are obtained by Blaizot and Krzywicki || and by Huang 
and Suzuki ||. Although the model is a simplified version of a more realistic situation, it 
captures the essential features of soft pion production such as the light cone singularity, the 
boost invariance and the longitudinal expansion. The electromagnetic spectra are found to 
be sensitive to the dynamical evolution of the classical pion field. Since the electromagnetic 
current coincides with the third component of the isovector current, the isospin angular 
oscillation of the condensate field can be a significant source for the electromagnetic emission 
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in the low mass (transverse momentum) region. 

The pionic part of electromagnetic interactions is introduced by replacing (we neglect 
the anomalous electromagnetic coupling) d^ir by d^iz + ^[tt, Q]A/j, in the kinetic part of 
lagrangian 

C em (x) = ie[7r~{x)d fl TT + {x) - Tx + (x)d^-(x)}A»(x) = -eJ^ l {x)A tM {x) , (1) 

where A^ is the photon field and J^ 1 is the third component of the classical isospin current 
Jp(x) = (7r(x) x <9 M 7r(a;))3. To 0(e 2 ), the matrix element for the dilepton emission from the 
classical field [pHTH|l is 

M = (£ + r\S^\0) = e 2 u(p 1 ) lu v{ P2 )D^(p 1 +p 2 )J^(p 1 +p 2 ) , (2) 

where J^{pi + p 2 ) is the Fourier transformation of J^(x), u and v are the Dirac spinors for 
the electron and positron, D Ufl is the photon propagator. The number of dileptons produced 
per final state phase space is 



dN e+e -=J2\M\ 



2 d?pi d 3 p 2 



(2tt)»2£ 1 (2vr)3 2 £ 2 
/ d*q5 4 (q - Pl - p 2 ) V \M\ 2 , ^ , rf3 f 2 • (3) 
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The dilepton differential distribution with respect to the lepton pair 4-momentum q is then 

d ~lFT = ^?[^-«V1^wrW, (4) 

where B — [1 + 2mf/g 2 ][l — Am 2 /q 2 ] 1 / 2 and is the lepton mass. Similarly, the direct 
photon momentum distribution is calculated 

dNry a 



;J cl »(q)J?(q) , (5) 



d?q/u q (2tt) 2 

where q stands for the photon 4-momentum. Moreover, splitting the current into parts 
parallel and orthogonal to q^\ J^(q) = q^J pa,ra {q) + J^(q), one finds that only the orthogonal 
part contributes to the emissions. The dependence of the electromagnetic spectra on the 
space-time history is entirely contained in the Fourier transformation of the isospin current. 



Although some more realistic numerical simulations on the classical pion field dynamics exist 
H and ideally one would like to perform the Fourier transformation on these solutions, we 
do not attempt to do it in this Letter. Instead, we examine some simple scaling solutions 
in the 1+1 boost invariant model which are first obtained by Blaizot and Krzywicki ||. 
Since the calculation is fully analytical, we hope to understand some qualitative features of 
dilepton production from the classical field. 
The lagrangian for the nonlinear a model is 



£ = iftr 
4 



f 2 r m 2 f 2 



d^{x)d^(x) + — ^tr tf(x) + , (6) 



where one defines the pion field by a + it • n = f^H with the constraint a = y f 2 — tz 2 . 
£ transforms like £ — > U^TiU^ under SU(2)l x SU(2)r rotations. We are interested in an 



idealized boost-invariant case |Tl] where the field 7r = n(r)7r(r, x±), i.e. the orientation n 
of the pion field is only function of the proper time r defined by r = \Jt 2 — z 2 and uniform 
in spatial rapidity and transverse coordinates x±. This would correspond to a single DCC 
domain in the transverse dimension. The transverse profile of the field amplitude tt 2 is 
determined by the shape and size of the DCC domain, which we shall take a Gaussian form 
as suggested in numerical simulations 0] 

g(x ± ) =exp[-\x ± \ 2 /R 2 D } , (7) 

where Rd is the size of a DCC domain (typically Rd = 1 ~ 2 fm 0). Such a class of 
solutions are well known to predict a distinctive neutral pion probability distribution 

dP 1 

(8) 



df 2v/J' 

where P is the probability of find a neutral pion fraction / = N^o/lN^o + N n ±). However, 
if there exists more than one DCC domain, as it is likely the case especially in heavy- 
ion collisions where the transverse dimension Ra of the interaction volume is much larger 
than the size of a DCC domain Rd, the field orientation n could depend on the transverse 
coordinates x±. In this case, the prediction (|j) in general no longer holds. Such a possibility 



can only be studied in detail in numerical simulations. In this Letter, we shall confine 
ourselves to the case when the single domain dominates and only comment on how our 
result might change in a multi-domain case. 

According to the "Baked Alaska" scenario [[[]] , the boundary conditions for the solutions 
are imposed such that the classical waves propagate forward in proper time, starting at an 
initial time To. In other words, we are interested in a retarded wave that vanishes when 
t < To and t < tq. Such a solution only exists when there is a source located on the 
hyperbola of equal t = tq. The "hot" source can be regarded as the summation of the 
degrees of freedom other than pion modes. Away from the source located near the surface 
of the light cone, the field propagates freely and one should have the (partial) conservation 
of (axial) isovector currents 

d„V» = d^TT x = , (r > r ) (9) 

d,A» = d^ad^iz - Trd^a) = -mlUiz . (r > r ) (10) 

For a function only of r, a partial derivative d^f{r) is equal to (x M /r)<i/ /dr where = 
(t, 2;, 0,0). Eq.(|9p yields the integration 

7T x d»n = fia^g{x ± ) , (11) 

where a is a dimensionless constant vector in isospin space. In the chiral limit, the axial 
current is also conserved, which defines another constant vector b in isospin space with a 
constraint a ■ b = || 

nd^a - adj* = flb^g{x x ) . (12) 

The pion field orientation vector n(r) precesses with proper time r around a in the plane 
perpendicular to a 0. 

The solution ( |TTD exhibits an inverse square singularity as r goes to zero. As pointed 
out in 0, the singular oscillation near r = does not result from the neglect of pion 
mass: treating the pion mass as a perturbation does not generate qualitative modification 
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of the solution at small value of r. Since the ratio x^/t can be regarded as the velocity 
of a comoving element with coordinate x^, the current diverges outward from the origin 
indicating the existence of a source near the space-time origin. The boundary condition is 
satisfied if one multiplies(|Tl"|) by a step function 9{r — tq). The isovector current in the whole 
space-time is 



V,(t, x ± ) = f 2 a-^9(r - r )g(x ± ) , (13) 



The source term is then calculated (note that d^x^fix)] = (d t , d z , 0, 0)[x ,J, f(x)]) 

d,V» = flaU(T-T )g( Xl _). (14) 

It is clear that the pionic part of isospin (also the electromagnetic) current is not conserved 
due to the existence of a source term at r = r while the integration constants a and b specify 
the orientation of the source in isospace or the initial condition. It is also straightforward 
to calculate the energy density associated with the internal chiral oscillations using the 
calculated classical currents V and A and the classical solution for the pion field [§| 

2/> 2 + b 2 ) 



T 2 



(15) 



As required by the chiral symmetry, the energy density depends only on the combination 
a 2 + b 2 . 

The dilepton and photon emissions probe the nature of the source using time-like mo- 
menta q = (go, Q\\, Qt) with q 2 = M 2 and q 2 = respectively. The Fourier transformation of 
the electromagnetic current can be readily calculated in the light cone variables 

j; l (Q) = fas J d 4 xe^9(r - r )g(x ± ) 

= -^/>3^ex P [-g^y4][J (M ± r ) -iN (M ± r )) , 



where g M = (go, q\\, 0, 0), M± = ^Jq 1 — gjj and Jo and iVo are the Bessel functions. For a 
direct photon, Mj_ = q?. In deriving (|T6|) , we have performed the r-integration from r to 



infinity. Since the electromagnetic current density decreases rapidly with the proper time, 
the integrated current turns out to be finite and the emission is most effective in the early 
evolution of the classical field. It is clear from fllCf ) that due to the singular behavior of 
the Neumann function for small argument, the emission spectra are sensitive to the initial 
time scale r . The exponential factor exp[— q^R 2 D /4] provides a strong q T cutoff so that the 
emission is only significant in low momentum region. 

The Fourier transformation of the electromagnetic current takes the form J^(q) = 
q ll f(M±, qx) where / is a scalar function, the orthogonal part is then 

J t ;(q) = (0,0,-q T )f(M ± ,q T ) , (17) 

and 

{ff(f-?Sr)tf{Q)J?{*) = q 2 T Ml\f(M ± ,q T )\ 2 (18) 

J^(q)JT(q) = -q 2 T \f(M ± ,q T )\ 2 . (19) 

The dilepton and photon distributions are calculated 

= g(^)^«^|jexp[-^y2][J o2 (M ± . ) + N!(M ± r o)] , (20) 
(/A " - %(7rRl) 2 fy 3 -ex P [-q 2 T Rl/2][J 2 (q T r ) + N 2 (q T r )} , (21) 



dydq T 8 v u ' "« 6 q T 

where q\ = M\ — M 2 for the dilepton. Both the dilepton and the photon spectra fall off 
exponentially for large transverse momentum which is the characteristic of the coherent 
production from a finite domain. The electromagnetic emission from a DCC domain is thus 
only important in the low momentum region where the spectra increase as some inverse 
power law of the momentum. 

The constant CI3 can be estimated using the symmetry argument. The important as- 
sumption in the disoriented chiral condensate scenario is the equal probability of all internal 
orientations. According to Eq. (|i~5|) , the energy density due to the internal chiral oscillations 
should be equally distributed among all components of vectors a and b. One thus has 
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^) = ^) = ¥B (< = 1,2,3). (22) 

where eo is the initial energy density carried by the internal chiral oscillationat at r = To, 
which may be roughly estimated assuming that it is comparable to the thermal energy 
density (in order for the DCC to have any physical significance). For some typical values 
e ~ 0.5T C 4 ~ 0.1 GeV/fm 3 and r = 5 fm/c, the average value of a\ is about 4.8. 
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Dilepton Invariant Mass M (MeV) 

FIG. 1. The dilepton invariant mass spectrum for different choices of initial time scale tq and 
the coherent field domain size Rr>. The initial energy density is assumed to be eo = 0.1 GeV/fm 3 . 
The typical thermal spectrum due to the 7r-7r annihilation is also plotted for a comparison 

In Fig.l we numerically integrate (|2~0"D over M± and plot the dilepton invariant mass 
spectrum for different values of To = 5, 10 fm/c and Rd = 1,2 fm. The 03 is calculated 
from (p^ ) assuming that the initial energy density for the internal chiral oscillations is 0.1 
GeV/fm 3 . Also plotted in Fig.l is the typical thermal spectrum of dilepton production 
mainly from tt-tt annihilations, which is taken from Ref. [|K|, where we have assumed the 
initial temperature of 160 MeV, the freeze-out temperature of 130 MeV, and the transverse 



8 



size of system Ra = 5 fm suitable for heavy-ion collisions. Most notably, there does not exist 
a pion mass threshold for the dilepton production from the coherent field: the spectrum rises 
even below M = 2m n . The finite pion mass plays no roles in the conserved isovector current 
whose third component directly couples to the photon. Therefore, the coherent field is most 
effective in producing lepton pairs in the low mass region M < 2m n . 
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Transverse Momentum q T (MeV) 



FIG. 2. The direct photon transverse momentum spectrum for different choices of tq and Rd- 
The initial energy density is assumed to be eo = 0.1 GeV/fm 3 . 

Plotted in Fig. 2 is the direct photon transverse momentum spectrum for different choices 
of r and Re- The outstanding feature in the spectrum is the sensitivity on the domain 
size parameter R D . The measurement of direct photon spectrum at very small transverse 
momentum can thus provide a sensitive probe of the spatial extension of the DCC domain. 

So far we have only considered a single DCC domain and its orientation oscillation in 
time. If the multi-domain configurations dominate the DCC production, the dilepton and 
photon can be coherently produced from different oriented regions. In particular, if there are 
two domains with opposite orientation (the molecule type), the singular rise of the spectra 
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near zero momentum will be suppressed and spectra at the moderate momentum will be 
enhanced. We plan to address these interesting questions in future publications. 

In conclusion, we have suggested that the study of electromagnetic emission can provide 
some information on the space-time history of a disoriented chiral condensate. The emis- 
sions due to the coherent isospin oscillation can be an important additional source to the 
electromagnetic signals from a dense hadronic matter, especially in low mass (momentum) 
region. It would be very interesting to see whether or not the qualitative features suggested 
in this paper can survive in more realistic models, and indeed may provide some mechanism 
for the observed low mass dilepton enhancement by the CERES and HELIOS collaborations 
at CERN SPS experiments [0. 

We would like to thank Hans Thomas Elze and Ina Sarcevic for useful discussions. This 
work was supported through the U.S. Department of Energy under Contracts Nos. DE- 
FG03-93ER40792 and DE-AC03-76SF00098. 
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